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Abstract: Supersymmetric (SUSY) models, even those described by relatively few pa¬ 
rameters, generically allow many possible SUSY particle (sparticle) mass hierarchies. As 
the sparticle mass hierarchy determines, to a great extent, the collider phenomenology of 
a model, the enumeration of these hierarchies is of the utmost importance. We therefore 
provide a readily generalizable procedure for determining the number of sparticle mass hier¬ 
archies in a given SUSY model. As an application, we analyze the gravity-mediated SUSY 
breaking scenario with various combinations of GUT-scale boundary conditions involving 
different levels of universality among the gaugino and scalar masses. For each of the eight 
considered models, we provide the complete list of forbidden hierarchies in a compact form. 
Our main result is that the complete (typically rather large) set of forbidden hierarchies 
among the eight sparticles considered in this analysis can be fully specified by just a few 
forbidden relations involving much smaller subsets of sparticles. 
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1 Introduction 


Low-energy supersymmetry (SUSY) [1] is a well-motivated paradigm for new physics be¬ 
yond the Standard Model, as it explains the lightness of the recently-discovered Higgs 
boson [2, 3], provides a dark matter candidate [4], and, in some models, allows for the uni¬ 
fication of gauge couplings, consistent with the low energy data [5]. That there is currently 
no compelling evidence of sparticle production at the CERN Large Hadron Collider (LHC) 
[6, 7] strongly motivates efforts to make sure that “no stone is left unturned” in searching 
for SUSY and/or setting more rigorous limits on the masses of sparticles. It is therefore 
important to understand all possible signatures of SUSY models; the rich phenomenology 
of SUSY means that there are unexplored “corners” of even well-known models like the 
CMSSM [8]. Such studies will become particularly valuable once a signal of new physics 
is observed, as they will help determine whether we are seeing a SUSY signal in the first 
place, and if so, what variety of SUSY model is responsible for it. 

However, the exploration of SUSY collider phenomenology in full generality is a very 
challenging endeavor, complicated by the multitude of scenarios and the large number of 
SUSY parameters. Several approaches have been tried: 

• Analysis of specific benchmark points [9-14]. The advantage of this approach is that 
it relies on well motivated and popular theory models. Benchmark points provide 
concrete, clear targets for discovery, and allow the sensitivities of different types of 
experiments to be compared in a meaningful way [15-17]. The downside is that 
the conclusions are very model-dependent and cannot be easily generalized to any 
arbitrary SUSY scenario. 

• Analysis of the (phenomenologically relevant) MSSM parameter space (pMSSM) [18- 
23]. The advantage here is that in principle one is exploring all corners of parameter 
space and encountering all phenomenologically interesting signatures. The disad¬ 
vantage is that the large dimensionality of the pMSSM parameter space makes full 
coverage impossible. This has prompted studies in a correspondingly smaller param¬ 
eter subspace, e.g., a nine parameter subset (pMSSM9) [24]. Alternatively, one could 
restrict one’s attention to a specific sector of the model, e.g., the four lightest states 
in the new physics particle spectrum [18, 25, 26]. (See also the simplified model 
approach below.) 

• Analysis of simplified models. The third, intermediate, approach is motivated by 
signature-based searches for new physics at the LHC. For any given experimental sig¬ 
nature, one may consider only the particles that are relevant to this specific channel, 
arriving at a so-called “simplified model” with only a handful of parameters [27, 28]. 
This approach inherits, to a degree, the advantages of the previous two methods, but 
the connection to the underlying high energy theory (and its fundamental parame¬ 
ters) is obscured. 

The difficulty in exploring the phenomenology of the pMSSM parameter space in full 
generality manifests itself in the “LHC inverse problem” [29] of mapping an experimentally 
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observed set of signatures at the LHC to a specific parameter point in theory space. The 
main challenge stems from the vastness of the SUSY parameter space — even under ideal 
circumstances the map is not unique and exhibits degeneracies. Furthermore, such maps 
so far have been constructed by (relatively sparse) scans in parameter space [30-33], and 
it is not clear to what extent the derived conclusions are robust and reliable. 

A possible resolution to these problems was put forward in Ref. [34], which proposed 
a more manageable and practical parameterization of theory space. The main idea is that 
the generic parameter space of n SUSY parameters can be thought of as the direct 
product of the set of all possible permutations Sn of those n parameters and the remaining 
coset MJ^/Sn- 

= Sn®R^/Sn. ( 1 . 1 ) 

The n parameters considered in Ref. [34] were all SUSY mass parameters, hence the per¬ 
mutations in Sn were named “mass hierarchies”, or “hierarchies” for short. The advantage 
of the decomposition (1.1) becomes evident when we consider the nature of a new physics 
signal at a collider and, in particular, the dependence on the two factors Sn and R^/Sn- In 
general, the collider phenomenology of any new physics model depends on certain qualita¬ 
tive and quantitative aspects, where the former are parameterized by Sn, while the latter 
depend mostly on R^/Sn- 

• Quantitative aspects. Those are the factors which determine the overall signal rate, 
such as: 

1. Signal production cross section. This depends on the magnitudes of the spar- 
ticle masses, which are encoded in the R^/Sn factor alone. Generally, heavier 
particles have smaller cross-sections and vice versa. 

2. Signal branching fractions. These are functions of the magnitudes of the mass 
splittings, which are also part of R^/Sn- For example, decay modes that are 
close to threshold, are kinematically suppressed, and the amount of suppression 
depends on exactly how close to threshold we are, i.e., on the size of the mass 
splitting. 

3. Signal acceptances and efficiencies. These depend on the hardness of the SM 
decay products observed in the detector, which in turn is again a function of the 
magnitudes of the mass splittings. Parameter space points in R"^/Sn with large 
mass splittings lead to harder leptons, jets, etc., and correspondingly higher effi¬ 
ciencies. On the other hand, smaller mass splittings lead to degenerate scenarios 
with lower efficiencies, where discovery becomes problematic, see e.g. [35-37]. 

4. Relative contribution of strong versus electroweak production. In the years lead¬ 
ing up to the LHC, it was usually assumed that strong SUSY production would 
dominate, and the first sign of SUSY would most likely be seen in squark and/or 
gluino production. Such expectations are based on the fact that strong produc¬ 
tion is enhanced due to the large strong coupling and the color multiplicity 
factors. However, this expectation is not true in models with heavy colored 
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superpartners and light electroweak superpartners, where strong production is 
kinematically suppressed relative to electroweak production. The amount of 
suppression depends on the relative size of the superpartner masses, which is 
again parameterized by W^/Sn- 

• Qualitative aspects. These are the factors that determine the type of discovery sig¬ 
nature one is looking for, namely, the identity and multiplicity of SM particles in the 
hnal state. These qualitative features are mostly determined by the mass hierarchy, 
and can be parameterized by an element of the Sn factor in (1.1). 

The crucial observation of Ref. [34] was that, as long as we are interested in the qual¬ 
itative aspects of SUSY collider phenomenology, we are justihed in focusing on Sn only 
and studying hierarchies of sparticles without any reference to the actual sizes of their 
masses. Note that the number of elements in Sn is a hnite number (namely, n!), which al¬ 
lows a fully exhaustive exploration and classihcation of the sets of experimental signatures 
associated with each hierarchy. Somewhat surprisingly, such studies were able to reveal 
previously overlooked corners of pMSSM parameter space with dramatic, yet relatively un¬ 
explored, multi-lepton signatures [34]. While the analysis of [34] was originally applied only 
to pMSSM9, it was subsequently extended to include third generation sfermions and/or 
i?-parity violation in the MSSM [38] and to the NMSSM [39]. 

In this paper we build on the work in [34, 38, 39] and demonstrate the cataloguing 
of ah possible hierarchies within specihc SUSY breaking models. Following [34], we shall 
consider only the sparticles shown in table 1. Their mass spectrum is parameterized by 
the parameters listed in the last row of the table. Some of these parameters correspond to 
several particles; for instance, the wino mass Mw describes the mass of both a mostly wino- 
like neutralino and a mostly wino-like chargino w^. (The remaining two gaugino mass 
parameters are denoted as Mb, the bino mass, and Mq, the gluino mass.) Likewise, we 
consider two degenerate generations of sfermions, so that each of the hve sfermion masses 
(Mq, Mb, Mo, Mb and Me) describes particles in both the hrst and second generation. 
Finally, in the case of left-handed squarks (sleptons) the parameter Mq {Mb) refers to both 
members of the isodoublet, ub and {cb and ub). Tree-level sparticle mixing and one- 
loop corrections [40] complicate this story somewhat, since the experimentally measured 
mass eigenvalues are slightly different from the soft mass parameters in table 1. In what 
follows, we shall assume that the soft mass parameters 

{Mq, Mu, Me, Mb, Me, Mb, Mw, Mg} (1.2) 

have been extracted from the data via a global ht along the lines of Refs. [41-44]. 

Given the eight experimentally measured mass parameters of table 1, the main question 
which we address in this paper is the following: can these measurements, by themselves, rule 
out specific SUSY-breaking scenarios? For our purposes here, a SUSY-breaking scenario is 
nothing more than a set of boundary conditions for the soft mass parameters of table 1 at 
some high energy scale. For concreteness, we shall illustrate our method with a SUGRA- 
inspired scenario [45-47] in which the boundary conditions are specihed at the GUT scale. 
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Table 1. The set of SUSY particles considered in this analysis, shorthand notation (symbolic 
and numeric) for each multiplet, and the corresponding soft SUSY breaking mass parameter. Note 
our slightly unconventional notation for the three gaugino masses: {Mg, M^y, Mq} instead of 
{Ml, M 2 , M 3 } — in this paper the latter are reserved for the corresponding boundary conditions 
(1.5-1.7) at the GUT scale. 


We shall consider several popular variations of this SUGRA model, in which one imposes 
different sets of boundary conditions: 

1. Scalar mass unification into SU{5) multiplets. Since the SU{5) is unbroken above 
the GUT scale MquT: we shall always assume that squarks and sleptons belonging 
to the same SU{5) multiplet have a common mass at Mcut'- 

Mq{Mgut) = Mu{Mgut) = Me{Mgut) = Mio, Ml{Mgut) = ^^{Mgut) = 

(1.3) 

In general, M 5 and Mio are independent input parameters — even if they start out 
equal at the Planck scale, they will be separated due to RGE running from the Planck 
scale down to the GUT scale [48, 49] or due to H-term contributions [50]. 

2 . Scalar mass unification into 50(10) multiplets. A more constrained version of the 
model arises if we assume 5O(10)-like unification, where 

M 5 = Mio = Mq. (1.4) 

3. Gaugino mass unification. In SU{5) GUTs, the gaugino masses may exhibit (some 
combination of) up to four^ different patterns [51, 52], so that in general 

Mb{Mgut) = Ml, (1.5) 

Mw{Mgut) = M2, ( 1 . 6 ) 

Mg{Mgut) = M 3 (1.7) 

can be taken as free parameters at the GUT scale. We shall also optionally consider 
the usual assumption of gaugino unification. 

Ml = A ^2 = M3 = Mij2. 

^One for each irreducible representation appearing in the symmetric product (24 x 24)s of two adjoints 
in SU{5). 
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Table 2. The GUT scale assumptions behind each of the eight different SUGRA models studied 
in the paper. For completeness, we also list the number, d, of GUT scale input mass parameters in 
each case, as well as the names of those parameters. 

4. Universal Higgs masses and third generation sfermions. Finally, although we are not 
explicitly considering the mass spectrum in the Higgs sector and the third genera¬ 
tion sfermions, those mass parameters feed into the renormalization group equations 
(RGBs) through the hypercharge trace 

3 

5 = Tr {YM^) = ^ (M^ - - Ml + Ml) , (1.9) 

i=l 

where the index i now runs over the three generations of sfermions. Extending (1.3) 
over the third generation as well, the last term in (1.9) is identically zero, and our 
last GUT scale assumption becomes 

S{Mgut) = 0 <;=> Mh^{Mgut) = Mh^{Mgut), ( 1 - 10 ) 

which is essentially the requirement of Higgs mass unification. 

In this paper, we will always assume the GUT scale unification of sfermion families in 
SU{5) multiples described by eq. (1.3). By either imposing or not imposing the remaining 
three GUT scale assumptions (1.4), (1.8) and (1.10), we obtain a total of 2^ = 8 different 
model scenarios which are listed in table 2. Each case has a different number of input mass 
parameters at the GUT scale (those parameters are listed explicitly in the last column of 
the table), so the sparticle spectrum is constrained to varying degrees in each scenario. 

We proceed to study the allowed hierarchies in each model, using the shorthand nota¬ 
tion from table 1 to label each hierarchy, ordering the particles from heaviest to lightest. 
Eor example, GQUDLWEB is a hierarchy with Mg > Mq > Mjj > Mjo > Ml > M\y > 
Me > Mb- In all eight models from table 2, the number d of input GUT-scale param¬ 
eters (listed in the fifth column) is less than the number of measured parameters (1.2). 
Therefore, the GUT-scale boundary conditions will impose certain relationships among 
the low-energy parameters (1.2), and as a result, depending on the specific model, some 
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hierarchies will be allowed, while others will be forbidden. The main purpose of this paper 
is to compile a complete catalogue of the allowed and forbidden hierarchies in each of the 
eight model scenarios from table 2. Our approach is completely general and can be easily 
applied to other SUSY-breaking mechanisms, where the soft mass parameters are given by 
a different set of boundary conditions at Mgut-, or are generated at an initial scale differ¬ 
ent from Mquj'. The method is not limited to SUSY models, and is equally applicable to 
non-SUSY scenarios, as long as the relevant mass parameters are evolved through a linear 
and homogeneous system of RGBs. 

This paper is one in a long line of works addressing the major question of how to test 
for supersymmetry once a signal of new physics is seen [53] . After the initial determination 
of (some of) the SUSY mass parameters, one would like to know whether the data is 
consistent with supersymmetry in general, or with a specific model (e.g., one of the models 
in table 2). Several approaches are possible: 

• Top-down approach. One can try to fit the data directly to the GUT-scale input 
parameters of the corresponding model [54-56]. In doing so, one is faced with the 
usual challenges of global minimization problems. In particular, one has to be careful 
to exhaustively cover all corners of parameter space, in order to be sure that a bad 
fit really excludes the model. 

• Bottom-up approach. Alternatively, one can use the measured SUSY parameters 
at the electroweak scale as boundary conditions and run the RGBs in reverse back 
to the GUT scale [57, 58]. This method provides a clear and intuitive picture of 
unification. However, since the RGB’s are coupled, one needs a sufficiently large 
number of measurements in order to completely specify the initial conditions at low 
energies. Therefore, at the initial stages of the discovery, when only a fraction of the 
SUSY mass spectrum has been measured, this method generally does not apply. 

• SUSY mass sum rules. The fact that the number of GUT scale inputs is less than the 
number of sparticle masses is reflected in the existence of certain model-dependent 
relations (“sum rules” for short) among the low energy parameters [59]. This idea 
has been richly explored in many different model scenarios [60-76], and has some 
relevance to our approach as well. Bor example, as illustrated below, we often find 
that the reason why certain mass hierarchies are not allowed is simply the fact that 
they violate one or more of the respective SUSY mass sum rules. However, our 
analysis will extend one step further and identify hierarchies which are consistent 
with the sum rules, but disallowed for other reasons — e.g., because they would 
require unphysical values for the GUT-scale boundary conditions or because they 
lead to tachyons in the spectrum. 

The main advantages of our approach compared to these earlier studies are the follow¬ 
ing: 

1. We often do not require knowledge of the complete mass spectrum in order to decide 
that a given model is ruled out. In the extreme cases, the knowledge of just two to 
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four mass parameters can be already sufficient to discredit a given model hypothe¬ 
sis, regardless of the values of the remaining mass parameters (which may even be 
unmeasured). 

2. We do not require any RGE analysis post-discovery, since all the required analytical 
work has already been performed ahead of time. 

3. Our results are robust and reliable in the sense that they are obtained analytically 
and do not rely on any scanning of parameter spaces or on numerical fitting. 

The paper is organized as follows. In section 2 we introduce our notations and de¬ 
scribe the usual numerical procedure which relates GUT-scale to weak-scale parameters. 
The next sections describe various ways of determining if and understanding why a hier¬ 
archy is allowed or forbidden. In section 3 we explore a linear algebraic way to identify 
allowed and forbidden hierarchies. In section 4 we present our results on the allowed sets 
of mass hierarchies within each theory model scenario from table 2. Then in section 5 
we demonstrate how the enumeration of forbidden hierarchies is enormously simplified by 
the consideration of forbidden sub-hierarchies. In sections 4.1 and 5.1 we examine the use 
of mass sum rules to understand intuitively the sets of allowed and forbidden hierarchies, 
employing the CMSSM as a concrete example. Section 6 is reserved for our conclusions. 
In appendix A we summarize our results on the sets of forbidden hierarchies. The Python 
code which contains all the hierarchy information and which accompanies this paper is 
described in appendix B. 

2 Notations and setup 

The connection between the phenomenological parameters (1.2) measured at the elec- 
troweak scale Mew and the corresponding input parameters (specified at the GUT scale 
Mgut) is provided by the MSSM RGEs, which allow exact analytical solutions at one loop. 
The starting point of our analysis is therefore the map [77] 



w= i (2,1) 


( 2 . 1 ) 



where IT is a vector of weak scale mass-squared parameters and is the corresponding 
vector of independent GUT scale mass-squared parameters for case m (see the last column 
in table 2). Note that the length of the vector and the meaning of its components 

both depend on the specihc SUSY scenario: 



dim(G(”"M (2.2) 
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where the number d of GUT scale input parameters is listed in the hfth column of table 2. 
The matrix , therefore, is an 8 x matrix which encodes the solution to the RGEs. 
It is instructive to illustrate eq. (2.1) with a couple of examples. 

First, consider the most general case of model m = 8, when there are = 6 inde¬ 
pendent inputs. The GUT scale parameter vector will then be^ 


Gf, Gf, Gf, Gf, Gf) = (m2o,m|,5,M|,M|,m2) , (2.3) 


where as usual we have rescaled the hypercharge trace parameter as 


S 


1 / «i 
66 \aG 



(2.4) 


Q' 

In what follows, the three a* = ^ parameterize the SM gauge couplings at the weak 
scale, while ac stands for the corresponding unified gauge coupling at the GUT scale. In 
principle, the precise values of the gauge couplings depend on the scale of sparticle masses. 
For TeV scale SUSY, however, we shall use the typical values ac — 0.041, ai ~ 0.017, 
02 — 0.034, and 0:3 ~ 0.118. 

For model m = 8 , then, Fq. (2.1) contains the 8 x 6 matrix 


1 U = 


( 


/ 

1 

0 

1 

C3 

C2 

36*^1 

\ 


m2 
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0 

-4 

C3 

0 

ici 



fMlo\ 

Ml 



0 

1 

2 

C3 

0 

Ici 



Ml 

Ml 



0 

1 

-3 

0 

C2 

hi 



S 

Ml 



1 

0 

6 

0 

0 

Cl 



Mi 

Ml 



0 

0 

0 

03 

0 

0 



Mi 

Mlv 



0 

0 

0 

0 

02 

0 



\MlJ 

\Ml) 

Mew 

1 

0 

0 

0 

0 

0 

Ol 

/ 



= ^( 8 ) (5(8) ^ (2.5) 


Mgut 


where the dimensionless coefficients a* and Cj are defined as 

a? 


CLi - n ^ 


a = 



(1 - ai). 


( 2 . 6 ) 


(2.7) 


Specifically, in our numerical analysis below we shall use the values oi = 0.170, 02 = 0.676 
and 03 = 8.29. Analyses performed after a SUSY discovery could use updated values and 
scan over the then current experimental uncertainties on these parameters. 

As another illustration, let us also discuss a simpler case — model m = 3 has three 
input GUT scale parameters: Mq, S and Mij 2 . A non-zero value for S at the GUT scale 

^Note that for the gaugino masses we continue to use the notation Mg, Mw and Mb- 
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can be generated by non-universal Higgs masses, thus in what follows we shall label this 
model as NUHM [78, 79]. The GUT scale parameter vector can be taken to be 


G(3) = Gf, Gf)) = [mI S, , 

and the corresponding 8x3 matrix 


1 

1 

C3 + C2 + ^Ci \ 

1 

-4 

C3 + gCl 

1 

2 

C3 + gCi 

1 

-3 

C2 + jCi 

1 

6 

Cl 

0 

0 

03 

0 

0 

02 

0 

0 

Ol / 


( 2 . 8 ) 


(2.9) 


Finally, let us also mention the case of CMSSM, in which there are only two GUT 
scale parameters: 

G(i) = (gS'\ G^')) = (m2, m2/2) , (2.10) 

and the corresponding 8x2 matrix given by 




/ 1 

C3 + C2 + ^Cl \ 

1 

C3 + |ci 

1 

C3 + gCi 

1 

C2 + iCi 

1 

Cl 

0 

03 

0 

02 

V 0 

Ol / 


( 2 . 11 ) 


3 Linear algebraic approach 
3.1 Theory 

We now develop a criterion to decide whether, within a given model m from table 2, a 
specific n-particle hierarchy is possible or not. Note that in general we consider values for 
n in the range 

2 < n < 8, (3.1) 

i.e., we study not only full hierarchies of all 8 particles in eq. (1.2), but also arbitrary 
subsets of less than 8 particles (n < 8). As an extreme example, which we shall use in this 
section to illustrate the formal math, consider the 2-particle hierarchy (re = 2) lUQ in the 
CMSSM model (m = 1). 

Recall that the weak-scale masses-squared were arranged in the vector W introduced 
in (2.1). Thus we can specify a given (sub-)hierarchy of interest by either a string of re 
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letters taken from the second row of table 1, or by a sequence of n integers (ii,i 2 , ■ ■ ■ ,in) 
denoting the corresponding components of the vector W and taken from the third row of 
table 1. For example, the hierarchy Mw > Mq can be represented by the letter string WQ 
as explained in the introduction, or by the integer pair (^ 1 ,^ 2 ) = (7,1) representing the 
respective components of the vector W, namely (Wj,Wi). Thus in general an n-particle 
hierarchy ■ ■ ■, in) is represented by an n-component vector 






(3.2) 


When we are dealing with a sub-hierarchy (n < 8), this vector has fewer than 8 components, 
and we find it convenient to always promote it back to an 8-component vector by adding 
the remaining W components Wj ^, Wj^ ■,■■■, '^js-n • 

H = , VFjj,, • • • , , Wji, lFj2, • • • ) ^i(8-n)) ’ ^ *2; ■ ■ ■ ) in}- (3-3) 


In other words, the first n components of H are the weak scale masses-squared in the 
n-particle hierarchy of interest (in decreasing order of mass), while the remaining 8 — 
n components are the remaining weak scale mass parameters, taken in arbitrary order. 
Operationally, we can build the vector H by multiplying W by the 8x8 square matrix 


%,1) 

hh,2) 

••• <^(*1,8) \ 

%,1) 

'^(*2,2) 

• • • <^(*2,8) 



• • • '^(in,8) 

<^01.1) 

^{j2,2) 


j{8 — n) 5 !) 

'^0'(8-n),2) 

^(i(8-n),8)/ 


(3.4) 


so that 


H = U^W. (3.5) 

For our example of the 2-particle hierarchy WQ = {Wt,Wi), we have ii = 7 and ^2 = 1. 
Then for concreteness we can take ji = 2, j 2 = 3, js = 4, jh = 5, = 6 and Jq = 8, so 

that eq. (3.5) reads 


(H,\ 


/o 0 0 0 0 0 1 o\ 
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W2 
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W3 

H3 


00010000 


1F5 


W4 

He 
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We 


W5 

Hj 


00000100 


Wr 


We 

\Hs) 


\oooooooy 


Ws/ 




(3.6) 


Viable mass spectra exhibiting the required hierarchy are those for which the first n 
components of the vector H are all positive and in decreasing order, while the remaining 8 — 
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n components are just positive (and in any order). Since we are interested in a hierarchical 
ranking of the hrst n components, it is convenient to define another vector, A„, in terms 
of the mass-squared differences: 


/ H1-H2 \ 

H2-Hs 


A„, = 


Hn 

H-n+l 


\ ^(8-n) / 


= DnH, 


where the 8x8 square matrix Dn is dehned by 


/I 

-1 

0 

... 0 

0 

0- 

■ 


0 

1 

-1 

... 0 

0 

0- 

• 0 


0 

0 

0 

... 1 

-1 

0- 

• 0 

_\6ij-6i+ij, for i < re; 

0 

0 

0 

... 0 

1 

0- 

• 0 

1 Sij, for i > n. 

0 

0 

0 

... 0 

0 

1- 

• 0 


\0 

0 

0 

... 0 

0 


1/ 



(3.7) 


(3.8) 


For the WQ example in the CMSSM, the corresponding matrix D 2 is 


(1 

-1 

0 

0 

0 

0 

0 

o\ 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

VO 

0 

0 

0 

0 

0 

0 

1/ 


Combining (2.1), (3.5) and (3.7), we get 

(3.10) 


A set of GUT scale parameters is able to produce a physically meaningful mass spec¬ 
trum for the re-particle hierarchy under consideration^ if and only if all eight components 
of A^™^ are positive. We have thus reduced our original problem of proving the existence 

n,H 

®Note that by construction depends on the hierarchy, so that the vector is defined on a case 

by case basis, as suggested by the index H. 
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of the hierarchy (3.3) to the problem of finding solutions such that all components of 
are positive. 

The mathematical problem, therefore, is to determine whether for the 8 x matrix 

= DnU n (3.11) 

n,H 


appearing in (3.10), there exists G 


such that the vector 


^ ( 3 _^ 2 ) 

n^H n,H ^ ' 

has only positive components. For this purpose, we can use a theorem known as Gordan’s 
Alternative, [80] which states that for every r x s matrix. A, either 


3 X G such that Ax ^ 0 


(3.13) 


or 

G M'" such that = 0, p > 0. (3.14) 

Here we have used the usual mathematical notation that for a vector a, a 0 indicates 
that all components of a are positive, while a > 0 only implies that all components are 
non-negative and in addition they cannot be all zero. 

We can now directly make use of Gordan’s alternative by identifying the matrix A with 
our matrix from (3.11). Then, if we can find a satisfactory solution to the equation 


71,H 

the corresponding hierarchy is forbidden, i.e., there are no values of the GUT scale pa¬ 
rameters such that ^ 0. It is often easier to determine whether a solution to 

n,H 

eq. (3.15) exists than to prove the existence of a vector for which 0. 

We shall demonstrate this with a few concrete examples below in Sec. 3.2. 

Note that this approach also allows us to add additional constraints on the range of 
the GUT scale parameters themselves. Since the GUT scale boundary conditions (the 
components of G*^™^) are masses squared, one could demand that they are positive as well. 
This is certainly true for the gaugino masses squared Mf, and optionally for the scalar 
masses squared M^g and M| [81]. The trace S can, of course, have either sign. 

If we demand that a particular GUT scale parameter G^™^ be positive, we simply 
supplement the matrix with a row vector e* with components {ei)j = 5ij. Then the 

matrix A becomes 


Then, 





n,H 


G 


(m) 


(3.16) 


(3.17) 


- 13 - 


Gordan’s Alternative tells us that either one can find acceptable^ values of the GUT scale 
inputs which will lead to the required hierarchy, or we can find a non-trivial solution to 
the matrix equation 


= 0, (3.18) 

such that all components of p are non-negative. Obviously, the procedure of enlarging the 
matrix A as in (3.16) can be repeated for as many GUT scale parameters as necessary. 

An alternative method, which may be easier to implement in practice, is to use the 
Findinstance [] command in Mathematica to obtain a list of all possible n-particle hier¬ 
archies in each of the eight SUSY scenarios. These lists are provided explicitly, along with 
methods to access this data, in the Python code described in appendix B. 


3.2 Specific examples 

3.2.1 A forbidden hierarchy: VUQ in the CMSSM (m = 1) 

Let us first apply the formalism of the preceding subsection to the case of the two-particle 
hierarchy WQ in the GMSSM model (m = 1). The matrix D 2 was given in (3.9), the 
matrix was dehned in (3.6), while was presented in (2.11). Thus from eq. (3.11) 
we can find the matrix 

2,H 


= 

2,H 


(3.19) 


/ — I 02 — C3 — C2 — •^Ci\ 

1 C 3 + C 2 + ^Ci 

1 C 3 + |ci 

1 C 3 -b |ci 

1 C2 + jCl 

1 Cl 

0 03 

Vo 01 / 

The GUT scale parameters are given in (2.10). Requiring both to be positive, the enlarged 
matrix A becomes 


A = 


/-I 
1 
1 
1 
1 
1 
0 
0 
1 
0 


a2 — C3 — C2 — -^CiX 
C3 + C2 + j^Ci 

C3 + |ci 
C3 + |ci 
C2 + jCi 
Cl 
03 
01 
0 
1 


/-I 
1 
1 
1 
1 
1 
0 
0 
1 
0 


-6.29\ 
6.97 
6.54 
6.49 
0.52 
0.15 
8.29 
0.17 

0 

1 / 


(3.20) 


'‘Meaning > 0. 
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Eqs. (3.18) now read 


Obviously, 



{pi,P2, ■ ■ ■ ,Plo) 


V 


1 

1 

1 

1 

1 

0 

0 

1 

0 


-6.29\ 

6.97 

6.54 

6.49 

0.52 

0.15 

8.29 

0.17 

0 

1 


= 0 . 


/ 1 \ 
0 
0 
0 
0 
0 
0 
0 
1 

y6.29y 


(3.21) 


(3.22) 


is a solution to (3.21) and therefore, by Gordan’s alternative, the hierarchy is not allowed. 
This can also be seen by noting that the first row of the A matrix (3.20) contains only 
negative components, and the first equation in (3.12) reads 


^ - 6.29 = -m2 - 6.29M2/2. 


(3.23) 


It is clear that the right-hand side cannot be positive, as long as both Mq and are 

positive as well. 


3.2.2 A forbidden hierarchy: LEU in the NUHM (m = 3) 

As a more complicated example of a forbidden hierarchy, consider the LEU sub-hierarchy 
in the CMSSM with non-universal Higgs masses [78, 79] (which in our notation simply 
means that S is no longer zero). For concreteness, we choose the remaining 5 particles to 
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be (5, D, G, W, and B and we form the extended vector H as in eq. (3.3): 


H = 


{ ^l\ 

Ml 

Ml 

Ml 

Ml 

Mg 

Mlv 

Wb) 


This defines the matrix 

/O 0 0 1 0 0 0 o\ 

00001000 
01000000 
^ _ 10000000 
00100000 
00000100 
00000010 
yO 0000001/ 

Using (2.9), (3.8) and (3.25) we find the corresponding matrix M' 1 to be 


/ 

0 

- 9 

C2 - fci^ 


/o 

-9 

0.37\ 


0 

10 

-C3 + fci 


0 

10 

-6.39 


1 

-4 

C3 + |ci 


1 

-4 

6.54 


1 

1 

C3 + C2 + T^Ci 

rs»/ 

1 

1 

6.97 


1 

2 

C3 + gCi 


1 

2 

6.49 


0 

0 

as 


0 

0 

8.29 


0 

0 

02 


0 

0 

0.68 

V 

0 

0 

ai/ 



0 

0.17/ 


(3.24) 


(3.25) 


(3.26) 


In this model scenario, we only demand that Mq and M^^ be positive, while S can have 
either sign. Therefore, the matrix (3.26) is extended with only two extra rows: 


/0 -9 0.37\ 

0 10 -6.39 
1 -4 6.54 

1 1 6.97 

1 2 6.49 

0 0 8.29 

0 0 0.68 

0 0 0.17 

1 0 0 

Vo 0 ij 


(3.27) 
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and the system of equations (3.18) to be studied is 


{pi,P2, ■ ■ ■ ,Pw) 



-9 

0.37\ 

0 

10 

-6.39 

1 

-4 

6.54 

1 

1 

6.97 

1 

2 

6.49 

0 

0 

8.29 

0 

0 

0.68 

0 

0 

0.17 

1 

0 

0 

VO 

0 

V 


= 0 . 


(3.28) 


This system admits strictly positive non-trivial solutions, for example 

P 2 = 0.9pi, pio = (0.9 c3 - C 2 + 0.25ci)pi ~ 5.38pi, 


(3.29) 


for any positive value of pi. According to Gordon’s alternative, the existence of the solution 
(3.29) implies that the hierarchy LEU is forbidden in the m = 3 model. We shall see below 
that LEU is indeed one of the 36 forbidden hierarchies in this model, see eq. (A.7). 

3.2.3 An allowed hierarchy: GBL in the NUHM (m = 3) 

We conclude this section with an example of an allowed hierarchy. Consider the sub¬ 
hierarchy GBL, which is one of the 198 allowed 3-particle hierarchies in model m = 3 (cf. 
table 3). As in the previous subsection, in model m = 3 the input parameters are still 
those given in (2.8), and the matrix is given by (2.9). For concreteness, we choose 
the remaining 5 particles to be Q, U, D, E, and W and form the extended vector H as in 
eq. (3.3): 

Ml 


which defines the matrix 


Uh = 


Mj 

m 2 
Ml 
Ml 

\Ml,J 


/o 0 0 0 0 1 0 o\ 
00000001 
00010000 
10000000 
01000000 
00100000 
00001000 
\0 0000010 / 


(3.30) 


(3.31) 
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Then, using (2.9), (3.8), (3.31) and (3.11), we obtain the enlarged matrix (3.16) 


/ 

0 

0 

as — ai^ 


/ 

0 

0 

8.12\ 


-1 

3 

ai — C2 — \ci 



-1 

3 

-0.35 


1 

-3 

C2 + jCi 



1 

-3 

0.52 


1 

1 

Cs + C2 + ^Ci 



1 

1 

6.97 


1 

-4 

cs + |ci 



1 

-4 

6.54 


1 

2 

Cs + |ci 



1 

2 

6.49 


1 

6 

Cl 



1 

6 

0.15 


0 

0 

02 



0 

0 

0.68 


1 

0 

0 



1 

0 

0 

V 

0 

0 

V 


1 

0 

0 

V 


The corresponding linear system of equations (3.18) is 


/ 0 
0 

(^ 8.12 


-1111110 
3 -3 1 -4 2 6 0 

- 0.35 0.52 6.97 6.54 6.49 0.15 0.68 



^ Pi^ 
P2 

\PioJ 


(3.32) 


(3.33) 


We want to show that this system has no solutions with p > 0. To easily see this, multiply 
the last equation by 4 and add to the first two equations, obtaining the equation 


32.47pi + 0.58 p2 + 0.09 p3 + 29.87p4 + 23.18p5 
+ 28.97p6 + 7.60 p7 + 2.70p8 + P9 + 4pio = 0. (3.34) 


As all of the coefficients in eq. (3.34) are positive, there is clearly no valid solution for p, 
and thus the hierarchy is allowed. 


4 Results: allowed n-particle (sub-)hierarchies 

For each of the eight model scenarios in table 2, only certain sets of hierarchies are allowed, 
in the sense that there exist values for the GUT scale parameters which will give rise 
to the given mass ordering at the weak scale. Note that in searching for viable hierarchies, 
we are interested not only in obtaining the masses in the specified ordering, but in addition 
we require that all masses are physical, i.e., that all mass squared parameters at the weak 
scale (2.1) are positive. This is why the hierarchy vector H was defined in eq. (3.3) in 
terms of all eight components of the vector W from (2.1). 

As indicated in (3.1), we shall also be interested in n-particle sub-hierarchies, i.e., sets 
of just a few particles with n < 8. In doing so, we are motivated by the experimental reality 
— during the phase of initial discovery of supersymmetry, it is very likely that only a few 
superpartners will be seen. Therefore, at that stage it would make no sense to ask questions 
involving the unseen yet superpartners. Instead, one should focus on the question, given 
what has been observed so far, is the data consistent with a given theory model or not? 
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Case 

n 

(m) 

2 

3 

4 

5 

6 

7 

8 

1 

36 

104 

190 

216 

148 

56 

9 

2 

42 

154 

382 

604 

570 

290 

61 

3 

49 

198 

519 

852 

827 

430 

92 

4 

51 

240 

757 

1536 

1887 

1252 

340 

5 

55 

312 

1277 

3232 

4560 

3254 

913 

6 

56 

330 

1521 

4806 

8684 

7824 

2699 

7 

56 

336 

1591 

5234 

10240 

10224 

3940 

8 

56 

336 

1648 

6028 

13778 

16502 

7766 

8!/(8-re)! 

56 

336 

1680 

6720 

20160 

40320 

40320 


Table 3. The number of allowed n-particle hierarchies for each of the eight models considered in 
table 2. 


We shall therefore allow n to vary within the full range (3.1). In general, for a given n, 
the number of all possible hierarchies which can be formed out of the set of eight observables 
(1.2) is given by and listed in the last row of table 3. The remaining rows of the 

table give the number of allowed n-particle hierarchies (in the sense described above), for 
each of the eight theoretical model scenarios from table 2. 

Table 2 reveals a significant reduction in the number of possible hierarchies. Consider, 
for example, the least constrained model. Case 8, where there are as many as 6 input 
parameters at the GUT scale. Even then, out of the 8! = 40, 380 8-particle permutations, 
only 7, 766 hierarchies are possible (a little over 19%). The other, more constrained models, 
exhibit a further reduction of the possible mass hierarchies. In the extreme case of the 
CMSSM (model m = 1 in table 2), there are only nine 8-particle possibilities, as explained 
in Section 4.1 below. 

When we consider sub-hierarchies with n < 8, the reduction in the number of allowed 
sub-hierarchies in table 3 is not as dramatic. This is because, when asking the question 
“Is this sub-hierarchy of n particles allowed or not?” we allow the remaining (unseen) 
8 — n particles to have arbitrary masses. In particular, in models m = 7 and m = 8, all 
possible 2-particle and 3-particle hierarchies are still represented in the table, in spite of 
the significant reduction in terms of the number of allowed 8-particle hierarchies. 

Of course, table 2 only tallies up the number of possible hierarchies in each model, but 
does not reveal which particular hierarchies are allowed. Given the large numbers seen in the 
table, giving here the complete list of all allowed hierarchies is impractical — the interested 
reader can easily generate the sets of allowed hierarchies from the accompanying code 
described in appendix B. This code can also answer queries about individual hierarchies 
(as well as sub-hierarchies) of interest. 

Nevertheless, one may still wonder whether it is possible to somehow display the com¬ 
plete information about the sets of allowed hierarchies which are hiding behind table 3. In 
the next section 5, we shall show that there is an elegant way of encoding and presenting 
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the same amount of information. The key idea is, instead of studying the number of al¬ 
lowed hierarchies (which is typically very large, as seen in table 3), to focus on the (much 
smaller) number of forbidden hierarchies, starting with the smallest possible values of n. 


4.1 The set of allowed hierarchies in the CMSSM 


Table 2 revealed that in each of the eight model scenarios considered, the number of GUT 
scale input parameters is less than the number (eight) of weak scale parameters under 
consideration. Hence the hierarchies must satisfy 8 — sum rules, which explains the 
reduction in the number of allowed hierarchies observed in table 3. 

We shall now illustrate how these sum rules help to identify allowed hierarchies using 
the example of the CMSSM (m = 1). In this case there are two input parameters (Mq and 


M‘^j2), so we have the following six sum rules: 

(5ci - 12 c 2 )(m 2 - m2) - 4ci(m2 - m2) = 0 , (4.1) 

3(5ci - 12 c 2)(M| - M^) - (7ci - 36c2 + 8603)(M^ - M^) = 0 , (4.2) 

3(5ci - 12 c 2)(M| - M^) + (20ci - 8603)(M^ - M^) = 0 , (4.3) 

(5ci - 12 c 2)M^ + 1203(M^ - M^) = 0 , (4.4) 

(5ci — 12 c 2)M^ -|- 1202 (Mq ~ ^u) ~ ^ > (^•^) 

(5ci - 12 c 2)M| -k 12oi(M^ - M^) = 0. (4.6) 


In obtaining these sum rules, we first used the first two rows of in (2.11) to express 
the GUT scale parameters. Mg and ™ terms of the weak scale parameters, Mq and 

Mh, 


Mn2 



4(4ci +^^ 3 ) _ Cl + 36 c2 + 36c3 ^^2 

15ci - 36 c2 ^ 15ci - 36 c2 ^ 


12 

5ci - 12 c2 


(m2 -m2). 


(4.7) 

(4.8) 


Then, the remaining six rows of the matrix result in the sum rules (4.1-4.6). 

We note that the overall mass scale does not have any impact on whether a given 
hierarchy is allowed or not. We can therefore remove one degree of freedom; here we 
choose this degree of freedom to be Mq. Since the CMSSM has only two input parameters 
(i.e., = 2), removing the overall scale as Mq leaves us with only one relevant degree 

of freedom, which we can take to be the ratio The sum rules (4.1-4.6) can then be 
rewritten as 


Mj) _ Cl - 12 c 2 / 4ci 
m 2 “ 5ci - 12 c2 Vm 2 y 5ci - 12 c2 ’ 

M| 8 ci — 36 c3 / M^j \ 7ci — 36c2 -k 8603 
m| “ 15ci -36 c2 \JR^) 15ci - 36 c2 


(4.9) 

(4.10) 
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M^/M^ 

E Q 

m"/m^ 

B Q 

-mW 

L Q 

M^/M^ 

W Q 

D Q 

M^/M^ 

U Q 

-M^/M^ 

Q Q 

M^/M^ 

G Q 


Figure 1. Dependence of the squared mass ratios (4.9-4.14) on the parameter We only show 
the physically viable range. 


Ml _ 35ci - 36 c2 - 36c 3 / \ 20ci - 36c3 

“ 15ci - 36c2 \J^) ~ 15ci - 36c2 ’ 

Ml ^ 12a3 f Ml \ _ 1203 

Mq 5ci — 12 c2 \Mq) 5ci — 12 c2 ’ 

_ 1202 (Ml\ 1202 

Ml ~ 5ci - 12 c2 \Ml) 5ci - 12c2 ’ 

Ml _ 12oi (Ml\ 12oi 

~Ml ~ 5ci - 12 c2 \Ml) ~ 5ci - 12c2 ’ 


(4.11) 

(4.12) 

(4.13) 

(4.14) 


which are clearly linear equations in In general, the ratio can take values in 

the interval (0,oo). However, not all values will lead to physically viable mass spectra, as 
illustrated in Figure 1. In particular, some of the weak scale masses squared can become 
negative and lead to tachyonic particles. As observed in in Figure 1, the most stringent 
such restriction arises from the constraint Ml > 0. Using eq. (4.11) this translates into 


Ml ^ 20ci — 36 c 3 
^ 35ci - 36 c2 - 36c3 


0.94. 


(4.15) 


In addition, we should also require positivity of the gaugino mass squared boundary con¬ 
dition (4.8), which in turn implies 


Ml 

Ml 


< 1 . 


(4.16) 


Requiring that Mg > 0 does not lead to any further restrictions, so that the physically 
allowed range for depicted in Figure 1, is (0.94,1). 
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For any given value of the hierarchy of weak scale masses is uniquely determined. 
Furthermore, whenever a pair of lines in figure 1 cross, the hierarchy changes, since the 
mass ordering of the corresponding particles is reversed. Therefore, we can enumerate the 
hierarchies by examining the values of the independent degree of freedom at which 
lines in Figure 1 cross. Let 


Ml 

Ml 


Mx=My 


(4.17) 


denote the value of at which the 


line and the 




line in Figure 1 intersect. It is 


easy to see that there are a total of eight intersection points, ordered as 


< 


Ml 

Ml 

V 

1 

II 

1 

< Mu 

Mw=Ml 

V 

1 

II 

i 

Ml 

< Mu 

< Mu 

< Mu 

Ml 

Mg=Mu 

Mg=Mb 

Mg=Ml 


< 


Mg=Mq 


Mg=Me 


(4.18) 


The eight points (4.18) divide the allowed interval into 9 sub-intervals, where each 
sub-interval corresponds to one specific hierarchy. This confirms that the CMSSM indeed 
has 9 possible 8-particle hierarchies, as shown in Table 3. Figure 1 reveals that they are 
(going from left to right): 


GQUDWLBE 

Mg 

> 

Mq 

> 

Mu 

> 

Md 

> 

Mw 


Ml 

> Mb 

> 

Me, 

(4.19) 

GQUDWLEB 

Mg 

> 

Mq 

> 

Mu 

> 

Md 

> 

M\y 


Ml 

> Me 

> 

Mb, 

(4.20) 

GQUDLWEB 

Mg 

> 

Mq 

> 

Mu 

> 

Md 

> 

Ml 

> 

M\y 

> Me 

> 

Mb, 

(4.21) 

GQUDLEWB 

Mg 

> 

Mq 

> 

Mu 

> 

Md 

> 

Ml 

> 

Me 

> M\y 

> 

Mb, 

(4.22) 

QGUDLEWB 

Mq 

> 

Mg 

> 

Mu 

> 

Md 

> 

Ml 

> 

Me 

> 

> 

Mb, 

(4.23) 

QUGDLEWB 

Mq 

> 

Mu 

> 

Mg 

> 

Md 

> 

Ml 

> 

Me 

> Mw 

> 

Mb, 

(4.24) 

QUDGLEWB 

Mq 

> 

Mu 

> 

Md 

> 

Mg 

> 

Ml 

> 

Me 

> M\y 

> 

Mb, 

(4.25) 

QUDLGEWB 

Mq 

> 

Mu 

> 

Md 

> 

Ml 

> 

Mg 

> 

Me 

> M\y 

> 

Mb, 

(4.26) 

QUDLEGWB 

Mq 

> 

Mu 

> 

Md 

> 

Ml 

> 

Me 

> 

Mg 

> Mw 

> 

Mb. 

(4.27) 


5 Results: the set of forbidden n-particle (sub-)hierarchies 


Table 3 may seem an unedifying melange of data. However, it turns out that if 
1 . one considers forbidden hierarchies rather than allowed hierarchies, and 


2. inspects a forbidden n-particle hierarchy for forbidden 2, 3, ..., n — 1-particle sub¬ 
hierarchies contained in it, 

then the picture becomes much simpler. Our results in terms of forbidden hierarchies 
are shown in table 4. The table is constructed as follows. For each model, we begin by 
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Case 

(m) 

n 

2 

3 

4 

5 

6 

7 

8 

1 

20 

0 

12 

0 

0 

0 

0 

2 

14 

0 

13 

0 

2 

0 

0 

3 

7 

36 

34 

10 

3 

0 

0 

4 

5 

18 

79 

14 

3 

0 

0 

5 

1 

6 

121 

168 

149 

12 

0 

6 

0 

6 

45 

216 

288 

98 

3 

7 

0 

0 

89 

176 

426 

434 

22 

8 

0 

0 

32 

148 

809 

398 

54 


Table 4. The number of irreducible forbidden n-particle hierarchies for each of the eight models 
considered in table 2. 


considering the 8 x 7 = 56 possible 2-particle hierarchies. For each of these 2-particle 
sub-hierarchies, we check whether it is present in at least one of the allowed 8-particle 
hierarchies counted in the last column of table 3.^ If the answer is “yes”, then the 2- 
particle sub-hierarchy is in principle allowed in this model, but if the answer is “no”, this 
2-particle sub-hierarchy will never appear in the model, and is therefore categorized as 
“forbidden”. The total number of forbidden 2-particle sub-hierarchies is then tallied up 
and listed in the n = 2 column of table 4, while the explicit lists of the forbidden 2-particle 
hierarchies appear in Appendix A, see eqs. (A.l), (A.3), (A.6), (A.11), and (A.13). Clearly, 
for any given model m, a specific 2-particle sub-hierarchy is either allowed or forbidden, 
so the total number ]^ 2 ^l,j.bidden forbidden n = 2 sub-hierarchies in the model and the 
total number of allowed n = 2 sub-hierarchies in the model add up to the total 

number of 2-particle sub-hierarchies; 


^^2Jorbidden ^ ^^2,allowed ~ 


(5.1) 


which is seen to hold for the results in the respective columns in tables 3 and 4. We also 
observe that as the number of GUT-scale inputs increases, the number of 

forbidden 2-particle sub-hierarchies decreases. In particular, for models 6, 7, and 8, we find 
that all 2-particle sub-hierarchies are possible. At the same time, in the most constrained 
model, the CMSSM, there are 20 forbidden 2-particle hierarchies, which will be discussed 
in detail in Sec. 5.1. 

Having thus determined the set of all forbidden 2-particle sub-hierarchies, we now move 
on to finding the set of all forbidden 3-particle sub-hierarchies. Clearly, a large number 
of 3-particle sub-hierarchies will be forbidden simply because they contain a forbidden 2- 
particle sub-hierarchy. Such cases are uninteresting to us, and we focus only on the search 
for irreducible forbidden 3-particle hierarchies, i.e., 3-particle sub-hierarchies which are 
composed of allowed 2-particle sub-hierarchies only. The number d^^^forbidden such “newly 

® Equivalently, we can directly check whether the 2-particle sub-hierarchy under consideration is one of 
the allowed 2-particle sub-hierarchies from the n = 2 column of table 3. 
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forbidden” 3-particle sub-hierarchies is listed in the third column of table 4. Interestingly, 
the results for ^^"flrbidden follow quite the same pattern that we already observed 

for ^ 2 ^lrbidden- particular, our 8 models fall into one of the following three categories: 

• Models in which in principle there exist forbidden 3-particle sub-hierarchies, but none 
of them is “newly forbidden”. In other words, all forbidden 3-particle sub-hierarchies 
already contain a forbidden 2-particle sub-hierarchy. According to table 4, models 
m = I and m = 2 are of this type. 

• Models in which there are no forbidden 3-particle sub-hierarchies to begin with. These 
are the models with a relatively large number of input parameters, namely models 
m = 7 and m = 8 . 

• Lastly, there are models in which we obtain non-trivial results at the level of 3-particle 
hierarchies. In models m = 3, m = 4, m = 5 and m = 6, there arise newly forbidden 
3-particle sub-hierarchies, which cannot be explained simply with the presence of a 
forbidden 2-particle sub-hierarchy. The number of “newly forbidden” 3-particle sub¬ 
hierarchies does follow the expected pattern — generally speaking, there are more 
such hierarchies for the models with fewer input parameters. 

Having accounted for the forbidden 2-particle sub-hierarchies and the newly forbidden 
3-particle sub-hierarchies, we then proceed to investigate the number ^^^flj-bidden newly 
forbidden 4-particle sub-hierarchies. The results are listed in the fourth column of table 4. 
Unlike the case of 3-particle sub-hierarchies, we hnd that now all models predict newly 
forbidden (4-particle) sub-hierarchies. In the case of models m = 7 and m = 8 , this is the 
first time we encounter forbidden sub-hierarchies. 

The process just described continues until we fill out the whole table 4. Comparing our 
results in table 3 for the allowed sub-hierarchies and the complementary results in table 4 
for the forbidden sub-hierarchies, we see significant simplihcation in the latter case in terms 
of counting hierarchies and bookkeeping. At the same time, it is worth emphasizing that 
tables 3 and 4 are simply two different ways to present our results on the restrictions 
imposed by various models on particle mass hierarchies. Whether we use the language of 
allowed or forbidden hierarchies, the amount of presented information remains the same. 

However, in practice there may be situations where one language is preferred over 
the other. For example, consider an initial discovery in which only a few supersymmetric 
particles have been observed. Let us focus on the simplest and most likely scenario where 
only two particles are detected: a lighter particle X with mass Mx, and a heavier partner 
Y with mass My {My > Mx)- Given just this single piece of information, what can one 
conclude about the various models from table 2? The answer is given in table 5, which 
lists all possible 2 -particle hierarchies, and the corresponding models in which a particular 
hierarchy is forbidden. There are several conclusions which can be drawn from the table: 

• Models m = 6, m = 7 and m = 8 are completely absent from table 5, which simply 
reflects the fact that in table 4 we did not find any forbidden 2-particle hierarchies 
for those models. This means that the discovery of just two supersymmetric particles 
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YX 

{My > Mx) 

X 

Q 

U 

D 

L 

E 

G 

W 

B 

Y 

Q 

- 








U 

1,2 

- 







D 

1 

1,5 

- 






L 

1 

1 

1,2 

- 





E 

1,2 

1,2 

1 

1 

- 




G 






- 



W 

1,2,3,4 

1,2 

1,2,3 



1,2,3,4 

- 


B 

1,2,3,4 

1,2 

1,2,3 

1,2 


1,2,3,4 

1,2,3,4 

- 


Table 5. The 56 possible 2-particle hierarchies YX (with My > Mx), where the heavier (lighter) 
particle is shown in the corresponding row (column) of the table. Each box in the table (except 
for the boxes along the diagonal) represents one of the 56 possible 2-particle hierarchies, and the 
entries inside each box indicate the model numbers m in which that hierarchy is forbidden. 

(no matter what they are) is insufficient to rule out SUSY models with large number 
of parameters like models m = 6, 7,8. 

• There are 36 2-particle sub-hierarchies which are present in all models — the 28 
cases above the diagonal in table 5, plus 8 more cases below the diagonal. If the 
data happens to be one of those 2-particle hierarchies, we will not be able to make 
any definitive statements about the eight theory models since they will all still be 
allowed. 

• The 2-particle hierarchies in which the heavier particle is colored while the lighter 
particle is not, are not restrictive at all. This is a well-known feature of SUGRA-type 
scenarios, in which the RGE evolution tends to split the spectrum in such a way 
that the colored superpartners are heavier than their electroweak counterparts. In 
particular, in all eight models, the gluino can be the heavier of the two particles. 
Similarly, in all models, the right-handed selectron can be the lighter of the two 
particles. 

• Among the most restrictive 2-particle hierarchies are those in which the lighter par¬ 
ticle is colored, while the heavier particle is not. For example, observing WQ, BQ, 
WG oi BG rules out four of the models right away. If both particles happen to be 
colored (or both are uncolored), then restrictive 2-particle hierarchies are those in 
which the lighter particle is charged under SU{2) while the heavier one is not, e.g., 
BW, BL, UQ, DQ. 

The python code described in appendix B not only generates all allowed hierarchies 
described in the previous section, but it can also easily reproduce the results on forbidden 
hierarchies discussed in the current section. For details on how to use the code and its func¬ 
tionality we refer the reader to appendix B. We shall conclude this section by considering 
a few illustrative examples which allow analytical treatment by means of sum rules. 
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5.1 Identifying all forbidden hierarchies in the CMSSM 

As seen in table 4, the CMSSM is characterized by 20 forbidden 2-particle sub-hierarchies 
and 12 newly forbidden 4-particle sub-hierarchies. These forbidden sub-hierarchies are 
explicitly listed in eqs. (A.l) and (A.2), respectively. Here we shall demonstrate how to 
re-derive those results with the help of the method of mass sum rules [59]. 

Our starting point is the relation between the weak-scale and GUT-scale parameters, 




/l 6.97 \ 

h4 


1 6.54 

Ml 


1 6.49 

Ml 


1 0.52 

Ml 


1 0.15 

Ml 


0 8.29 

^<w 


0 0.68 

\kj 


\0 0.17/ 



(5.2) 


The last three equations express the gaugino masses in terms of a single parameter, M^j^- 
Noting the different proportionality coefficients, it is easy to see that the CMSSM predicts 
the gaugino mass hierarchy 


Mq > M'^ > 


(5.3) 


This observation immediately rules out the three 2-particle hierarchies 


{BW,WG,BG}. (5.4) 


Another way to arrive at this result is to note that one can eliminate the M 1/2 parameter 
and obtain two independent gaugino mass sum rules, e.g. 


0.68M^ 

- 8.29M^ 

= 0, 


(5.5) 

0.17M^ 

- 0.68M| 

= 0, 


(5.6) 

which can be equivalently written as 






- 0.68(M^ - 

m2) 

- 7.61M^ 

= 0, 

(5.7) 

-0.17(M| - 

M^) 

- 0.51M| 

= 0. 

(5.8) 


Note that the numerical coefficients on the left-hand side of those equations are all negative. 
Thus the hierarchy WG, where > Mq, violates the hrst sum rule, eq. (5.7), while the 
hierarchy BW, where Mb > My/, is in contradiction with the second sum rule, eq. (5.8).® 
Similar logic can be applied to the sfermion sector, where the dependence on the 
parameter Mg is the same for all sfermion masses squared, and the mass splittings are 

®The impossibility of the third forbidden gaugino mass hierarchy, BG, follows from a linear combination 
of eqs. (5.5) and (5.6), e.g. 

0.17Mg - 8.29M| = -0.17(M| - Ma) - 8.12M| = 0. 
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induced only due to the different dependence on The CMSSM predicts the scalar 

masses in the order 

m2 >Ml>Ml>Ml> Ml, (5.9) 

which rules out ten additional 2-particle hierarchies, namely 

{EL, ED, EU, EQ, LD, LU, LQ, DU, DQ, UQ} . (5.10) 

This conclusion can also be justified with suitable sum rules, e.g. EL is ruled out by the 
relation 

-0.17(M|-M|)-0.37M| = 0. (5.11) 

Eqs. (5.4) and (5.10) account for 13 of the 20 forbidden 2-particle hierarchies. The 
remaining seven forbidden hierarchies, namely 

{WQ,BQ,WU,BU,WD,BD,BL} (5.12) 

can also be understood in terms of sum rules involving two sfermion masses and one gaugino 
mass: 


WQ : 
BQ : 
WU : 
BU : 
WD : 
BD : 
BL : 


-6.82(M^ - M^) - 6.14M^ - 0.68M| = 0, (5.13) 

-6.82(M| - M^) - 6.65M^ - 0.17M| = 0, (5.14) 

-6.39(M^ - Ml) - 5.71MI - 0.68M| = 0, (5.15) 

-6.39(M| - Ml) - Q.22MI - 0.17M| = 0, (5.16) 

-6.34(M^ - Ml) - 5.66MI - 0.68M| = 0, (5.17) 

-6.34(M| - Ml) - 6.17MI - 0.17M| = 0, (5.18) 

-0.37{mI - Ml) - 0.20MI - 0.17M| = 0. (5.19) 


As indicated in table 4, in the CMSSM model the forbidden 2-particle hierarchies, 
eqs. (5.4), (5.10) and (5.12), completely specify the forbidden 3-particle hierarchies as 
well. However, there still remain twelve newly forbidden 4-particle hierarchies which do 
not contain any forbidden 2-particle hierarchies and need to be motivated by a different 
argument. 

The idea is the following. Recall from Fig. 1 that in the allowed range for -r^ there 
exist eight crossing points (4.17) where we switch from a hierarchy XY (to the left of the 
crossing point) to a hierarchy YX (to the right of the crossing point). Now consider two 
such crossing points, XiYi —)■ YiXi and X 2 Y 2 —)• Y 2 X 2 , where all four particles Xi, X 2 , Yi 
and Y 2 are different. For definiteness let us also assume that 


M| 

K 


< 




K 


Mx^^My^ 


Now consider the 4-particle hierarchy 12 ^ 2 ^ 1 ^!. The XiYi bit requires 


Mu 

Ml 


Ml 

Ml 


Mx^=My^ 


(5.20) 


(5.21) 
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while the Y 2 X 2 bit in turn requires 


Mi 

Ml 


Mx2 =^Y2 


K 


(5.22) 


In light of our original assumption (5.20), eqs. (5.21) and (5.22) place contradictory re- 
quirements on the parameter -rX, thus the hierarchy Y 2 X 2 X 1 Y 1 is not allowed. 

From (4.18) we find a total of 13 pairs of crossing points for which all four particles 
involved are different. This gives us 13 possible candidates for newly forbidden 4-particle 
hierarchies: 


{QGBE, QGWL, QGWE, UGBE, UGWL, UGWE, 

DGBE, DGWL, DGWE, LGBE, LGWE, LWBE, EGLW} . (5.23) 

However, the last 4-particle hierarchy, EGLW, contains a forbidden 2-particle sub-hierarchy, 
EL, and is therefore not newly forbidden. This leaves us with exactly 12 newly forbidden 
4-particle hierarchies, which are precisely those listed in (A.2). 

5.2 Examples of forbidden hierarchies in the most general case (m = 8) 

In the previous subsection (5.1) we discussed the classification of forbidden hierarchies in a 
simple case like the CMSSM model. We shall now consider the other extreme, namely, the 
most general case of model m = 8, where all 6 GUT scale input parameters are a priori 
unconstrained. According to table 4, in the case of model m = 8, there are 32 forbidden 
4-particle hierarchies, which are listed in eq. (A. 18). In this subsection we shall illustrate 
a few specific examples. 

As we saw in the case of the CMSSM, a mass hierarchy may be forbidden because it 
violates a mass sum rule. In the case of model m = 8, there are 6 input model parameters, 
which enforces two mass sum rules among the 8 observable mass parameters, e.g. 

M‘^ +Ml-Ml- Ml- 1.56M^ 0.99M| = 0, (5.24) 

Ml + 2Ml - 2Ml -Ml- 2MMl + 1.44M^ 0.74MI = 0 , (5.25) 

5.2.1 The mass hierarchy GWQB 

Let us hrst consider the forbidden 4-particle hierarchy GWQB and rewrite the sum rules 
(5.24) and (5.25) in terms of the positive mass differences, e.g. Mq — M^, — Mq and 

Ml-Ml: 

- 1.56(M^ - Ml) - 0.56(M^ - M|) 0.42M| + mI - mI - mI = 0, (5.26) 

-2.34(M^ - M^r) - 0.91(M^ - M|) - 0.17M| + M^ + 2MI - 2MI - MI = 0. (5.27) 

If this manipulation renders one of the sum rules (or a linear combination of them) in a form 

where all terms have numerical coefficients of the same sign, the hierarchy will be clearly 
forbidden, due to the sum rule. Unfortunately, this is not the case here — the coefficients 
in both (5.26) and (5.27) have alternating signs, and the mass parameters can be suitably 


- 28 - 



adjusted to make the left-hand sides of those equations vanish. Therefore, the sub-hierarchy 
GWQB is not forbidden by the sum rules alone, and we need to further investigate the 
GUT scale boundary conditions G®. For this purpose, we shall invert eq. (3.12) and solve 
for G® in terms of the weak-scale mass parameters. Because of the two sum rules (5.24) 
and (5.25), only 6 of the weak-scale mass-squared parameters are linearly independent, 
and one possible choice would be to supplement the masses of the 4 particles entering the 
hierarchy GWQB under consideration, with the two slepton masses squared and MJ,. 
Then the relevant = 6 equations from the system (3.12) are 


(Ml - M^\ 


/I -1 0 0 0 0\ 


/O 0 0 0 0 1 0 0\ 





0 1-10 0 0 


00000010 


Ml 

Ml-Ml 


0 0 1-10 0 


10000000 


s 

Ml 


00 0 100 


00000001 


Mi 

Ml 


0 0 0 0 1 0 


00010000 


Ml 

\ Ml ) 


\0 0 0 001/ 


\0 0 0 0 1 0 0 0/ 


\Ml ) 


(5.28) 


Inverting this equation, we find 


fMlA 


/-0.94 

-1.8 

-0.6 

-0.45 0 

-0.2 \ 

Ml 


0.47 

0.18 

-0.42 

-1.16 1 

0.6 

s 


0.16 

0.3 

0.1 

-0.07 0 

0.2 

Mi 


0.12 

0.12 

0.12 

0.12 0 

0 

Ml 


0 

1.48 

1.48 

1.48 0 

0 

\MlJ 


V 0 

0 

0 

5.89 0 

0 / 


(Ml - M^\ 

Ml-Ml 

Ml 


\ 


M? 

Mi 


/ 


(5.29) 


where we have substituted numeric values. Note, in particular, the first equation, which 
specifies the value of the GUT scale parameter M^q as 

Mfo = -0.94(M^ - M^) - 1.8(M^ - M^) - 0.6(M^ - M|) - 0.45M| - 0.2M|. (5.30) 

Since all mass terms on the right-hand side have negative coefficients, the weak-scale masses 
squared are positive by definition, and the mass-squared differences are positive by the 
assumption of the hierarchy GWQB, we conclude that the GWQB hierarchy necessarily 
requires a tachyonic value for M^q and is therefore forbidden. 

Recall that above we had the freedom of choosing two additional weak-scale parameters 
to form the vector in the left-hand side of (5.28). It is instructive to see what happens if 
we had chosen a different set of weak-scale mass parameters, say and M£ instead. In 
that case eq. (5.28) is replaced with 


(Ml - M^\ 


(1 -1 0 0 0 0\ 


/O 0 0 0 0 1 0 0\ 


fMl,\ 

Mw - Mq 


0 1-1 0 0 0 


00000010 


Ml 

Ml - Ml 


0 0 1-10 0 


10000000 


s 

Ml 


0 0 0 100 


00000001 


Mi 

Ml 


0 0 0 0 1 0 


01000000 


Ml 


0 0 0 0 1 / 


VMf/ 


(5.31) 
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and inverting, we get 


/Mlo\ 


/-0.78 -1.36 -0.56 -0.66 0.2 0\ 


/Ml - Mlr\ 

Mi 


0 -1.15 -0.55 -0.55 -0.6 1 


Mw ~ Mq 

S 


0 -0.14 0.06 0.13 -0.2 0 


Ml-Ml 

Mi 


0.12 0.12 0.12 0.12 0 0 


Ml 

Mi 


0 1.48 1.48 1.48 0 0 


Ml 

\Ml ) 


\ 0 0 0 5.89 0 0/ 


\ Ml J 


Note that the matrix on the right hand side does not have all negative entries in any given 
row, so the contradiction is not immediately obvious. However, recall the existence of the 
two sum rules (5.26) and (5.27), which define the remaining two parameters and Mj,: 


Ml 

Ml 


( 0.78 -0.66 0.34 0.59 -1 1 

y - 0.78 - 2.22 - 0.22 1.02 -1 0 


(Ml-M^\ 
Ml-Ml 


Ml 

Ml 

\ Ml J 


Putting together (5.32) and (5.33), we obtain 


[Mlo\ 


/-0.78 -1.36 -0.56 -0.66 0.2 0\ 


Mi 


0 -1.15 -0.55 -0.55 -0.6 1 


/Ml - M^\ 

S 


0 -0.14 0.06 0.13 -0.2 0 


Mw ~ Mq 

Mi 


0.12 0.12 0.12 0.12 0 0 


Ml-Ml 

Mi 


0 1.48 1.48 1.48 0 0 


Ml 

Ml 


0 0 0 5.89 0 0 


Ml 

Ml 


0.78 -0.66 0.34 0.59 -1 1 


\ Ml j 

\Ml) 


\-0.78 -2.22 -0.22 1.02 -1 Q ) 



From a linear combination of the first and last equation, we obtain; 


(5.33) 


(5.34) 


2MIq + mI = -2.34(M^ - Mir) - 4.93(M^ - M^) 
- 1.33(M^ - Ml) - 0.3MI - 0.6M^. 


(5.35) 


This sum rule cannot be satisfied because the LHS must be positive while the RHS is clearly 
negative. The lesson is that, depending on our choice of parameters for the inversion, we 
may have to use the sum rules (5.26) and (5.27) in order to show that a hierarchy is 
forbidden. 
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5.2.2 The mass hierarchy WQGL 

As another example, consider the 4-particle sub-hierarchy WQGL. The RGE solutions are 




/-I 0 -1 -6.48 0.19 -0.004\ 


(Mix 

m2 -m2 


1 0 1 -1.81 0.49 0.004 


Ml 

m2 - Ml 


0 -1 3 8.29 -0.49 -0.04 


S 

Ml 


0 1 -3 0 0.49 0.04 


Mi 

Ml 
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Mi 

\ Ml ) 


\ 1 0 6 0 0 0.15/ 


\MlJ 


Selecting and as additional parameters to invert the RGE equation, we have 


/MlA 


/-0.69 0.27 -0.57 -0.57 0.12 -0.08 \ 


fM^r Mq\ 

Ml 


-1.15 -0.55 -0.55 0.45 -0.6 0 


m2 -m2 
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Ml 
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Ml 

\Ml ) 


\ 6.86 -2.68 1.05 1.05 4.77 4.77/ 


\ Ml j 


Again, we do not have a single row with negative coefficients, but forming the sum of M^q 
and we find 


Mfo + Mi= -1.84 (M^-m2)-0.28 (m2 -m2) 

-1.12(M^-M|) - 0.12 M|-0.48 M^-0.08 M|. (5.38) 


This sum rule cannot be satisfied, since all terms in the RHS have negative coefficients. 
Therefore, the hierarchy WQGL is forbidden. 


6 Conclusions and summary 

The search for SUSY is the paramount experimental challenge for Run II at the LHG. 
Once SUSY is discovered, it may provide invaluable clues about GUT-scale physics. The 
measured pattern of SUSY particle masses will play an important role in this quest, as we 
have explained above. Our proposal is to consider the relative ordering, or the “hierarchy”, 
of the measured sparticle masses. By analyzing the one loop SUSY RGE’s, it is relatively 
straightforward to derive the complete set of allowed hierarchies for a given choice of GUT 
scale boundary conditions. In this paper we considered hierarchies involving up to eight 
weak scale masses (first/second generation masses for each of the five sfermion families, and 
three gaugino masses), and analyzed eight different GUT-scale model scenarios (table 2). 
Our results are listed in appendix A and can be reproduced with the accompanying Python 
code described in appendix B. We also provided some intuitive arguments, based on mass 
sum rules and linear algebra tricks, to better understand and justify those results. 

The advantage of the approach presented in this paper is that it allows to draw defini¬ 
tive conclusions based on only partial information. In particular, we have seen that the 
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Figure 2. A schematic view on how hierarchies shed light on GUT-scale boundary conditions. 
Once a set of SUSY particle masses is observed (left), GUT-scale boundary conditions (center) 
incapable of producing this sparticle mass hierarchy are ruled out by querying our database. This 
conclusion can be independently verified by suitable mass sum rules or explicit running of the RGBs 
(right). 

knowledge of a sub-hierarchy of a few sparticle masses, measured in the very early days 
of discovery, is sometimes sufficient to rule out a specific model. Consider, for example, 
the most general model, m = 8, which has 6 GUT-scale input parameters. In principle, 
one would need 6 independent measurements in order to fully reconstruct the GUT-scale 
physics. On the other hand, we have shown that in this model, there are as many as 32 
forbidden 4-particle sub-hierarchies. Therefore, the m = 8 model can be ruled out with 
only 4 suitable measurements, if the data happens to point to one of the 32 forbidden 
hierarchies. This procedure is pictorially illustrated in figure 2. 

The analysis presented here can be extended in many directions. Obviously, one may 
consider more of the MSSM particles, or generalize to non-minimal models like the NMSSM. 
Our approach is not limited to SUGRA-type boundary conditions and can be applied in 
other SUSY breaking scenarios as well. Armed with the toolkit of allowed and forbidden 
hierarchies, we look forward to analyzing GUT scale models after a SUSY discovery at the 
LHC! 
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A Results on forbidden hierarchies 


In this appendix, we provide explicit lists of “irreducibly-forbidden” or “newly-forbidden” 
re-particle sub-hierarchies in various GUT scale scenarios (“models”)'. These are the for¬ 
bidden re-particle hierarchies which do not contain a forbidden sub-hierarchy. 

A.l Model m = 1 (CMSSM) 

The forbidden (and hence irreducibly-forbidden) 2-particle hierarchies are 

WU, BW, DQ, BL, ED, BU, WQ, LQ, EQ, EL, 

LD, DU, WD, BG, WG, BQ, LU, EU, BD, UQ. (A.l) 

The irreducibly-forbidden 4-particle hierarchies are 

UGBE, QGWE, UGWL, LWBE, QGWL, DGBE, 

DGWL, QGBE, LGWE, LGBE, UGWE, DGWE. (A.2) 

The are no newly forbidden re-particle hierarchies for re > 5 in this model. It is astonishing 
that all the information about the 8! possible hierarchies in the CMSSM is encoded in the 
two lists (A.l) and (A.l). 

A.2 Model m = 2 

For model 2, the forbidden 2-particle hierarchies are 

WU, BW, BL, BU, WQ, EQ, LD, WD, BG, WG, BQ, EU, BD, UQ. (A.3) 

The newly forbidden 4-particle hierarchies are 

GDLU, UGBE, DQWL, QGWE, GDLQ, GUED, DUEL, 

DGWL, QGBE, UGWE, QDLU, DQEL, GQED-, (A.4) 

and the newly forbidden 6-particle hierarchies are 

UDLWBE,QGDUWL. (A.5) 

There are no other newly forbidden hierarchies. 

A.3 Model rei = 3 (NUHM) 

For this model, the forbidden 2-particle hierarchies are 

BW,WQ,WD,BG,WG,BQ,BD-, (A.6) 

the newly-forbidden 3-particle hierarchies are 

DUQ, UEQ, EQD, LUQ, UDQ, LEU, GDQ, EUD, ELQ, LUD, ELD, LQU, 

^For the definition of the eight models, see table 2. 
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LEQ, EUQ, WEU, DGQ, LED, LUE, BEU, LGQ, DLQ, WUE, EGQ, LGD, 
EGD, QED, LDQ, BUE, GLD, GED, UED, LDU, BLE, GLQ, BEL, GEQ; (A.7) 

the newly forbidden 4-particle hierarchies are 

DQBE, GDLU, UGBL, DEBU, UQWL, DUBE, QDEU, DELU, GQLU, 

GDWU, GQBU, GDEU, UGWL, DEWU, QEWU, WEBL, GEBU, QELU, 
GQEU, GDBU, GQWU, GELU, QEBU, UDBL, QDWU, UDEL, UQEL, 

QDLU, GEWU, QULD, DQWE, QDBU, UQBL, DQLE- (A.8) 

the newly forbidden 5-particle hierarchies are 

QGWEL, QGUWE, DGWEL, QGUBE, QGWLE, 

DGUWE, QUGWE, DUGWE, DGWLE, QUGBE-, (A.9) 

and the newly forbidden 6-particle hierarchies are 

QGUDWL,QDEGWL,QDEGBL. (A.IO) 

There are no newly forbidden 7-particle or 8-particle hierarchies in this model. 

A.4 Model m = 4 

For this model, the forbidden 2-particle hierarchies are 

BW,WQ,BG,WG,BQ-, (A.ll) 

and the newly-forbidden 3-particle hierarchies are 

UEQ, WDU, BLD, BUD, EUQ, WEU, BEU, WUE, EGQ, 

BUE, WLD, BDU, BLE, WUD, BEL, BDL, WDL, GEQ. (A.12) 

Additionally, as shown in table 4, there are 79 newly forbidden 4-particle hierarchies, 14 
newly forbidden 5-particle hierarchies, and 3 newly forbidden 6-particle hierarchies. (There 
are no newly forbidden 7-particle or 8-particle hierarchies.) In order to save space, we do 
not list these hierarchies explicitly here, however they can be easily obtained from the 
python code described in appendix B. 

A.5 Model m = 5 

There is only one forbidden 2-particle hierarchy 

DU, (A.13) 

and six newly-forbidden 3-particle hierarchies: 

BGU,GWQ,WGQ,EGU,GBU,GEU. (A.14) 

Additionally, as shown in table 4, there are 121 newly forbidden 4-particle hierarchies, 168 
newly forbidden 5-particle hierarchies, 149 newly forbidden 6-particle hierarchies, and 12 
newly forbidden 7-particle hierarchies. (There are no newly forbidden 8-particle hierar¬ 
chies.) In order to save space, we do not list these hierarchies explicitly here, however they 
can be easily obtained from the python code described in appendix B. 
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A.6 Model m = 6 


In this model, all 2-particle hierarchies are allowed. The forbidden 3-particle hierarchies 
are: 


BGU,GWQ,EGU,GBU,GEU,WGQ. (A.15) 


The newly-forbidden 4-particle hierarchies are 


UEWQ, WEUQ, WDLQ, LQGD, WBUQ, LGQD, UQWE, WQDL, WULQ, 
UQWL, DQWL, UWQE, GDLQ, WQUL, UWEQ, GUED, WUQL, UBWQ, 
WDQL, UWQB, GUBQ, LGDQ, GQWU, WUEQ, GQWD, BUGE, UQWB, 
WUQB, DLWQ, DWQL, DWLQ, GLQD, EUGB, WUQE, UWBQ, UWLQ, 
WUBQ, WLDQ, UWQL, GUEQ, GUBD, WQGL, GQLD, GLDQ, WLUQ. (A.16) 


Additionally, as shown in table 4, there are 216 newly forbidden 5-particle hierarchies, 
288 newly forbidden 6-particle hierarchies, 98 newly forbidden 7-particle hierarchies, and 
3 newly forbidden 8-particle hierarchies. In order to save space, we do not list these 
hierarchies explicitly here, however they can be easily obtained from the python code 
described in appendix B. 


A.7 Model m = 7 

In model 7, all 2-particle and 3-particle hierarchies are allowed. The forbidden 4-particle 
hierarchies are: 

LQGU, WDQU, BGUE, GWQB, GULQ, WDEQ, DLQU, DWQB, DULQ, WDLQ, 
LUGQ, DBEU, WGQU, DEBU, GLWQ, WULQ, LGUQ, UQWL, WDBQ, DWUQ, 
WGLQ, DUBE, DBWQ, WQUL, GQLU, BEDU, WUGQ, LDEQ, GBUE, BGEU, 
BDUE, WGQB, GBWQ, WUQL, BEGU, WDQL, GEBU, GBEU, WDQB, GWBQ, 
GQWU, LUDQ, DEWQ, GUWQ, DLEQ, DELQ, DWQE, DWEQ, GWUQ, BUGE, 
GWQL, WEDQ, WUDQ, DLWQ, WGQL, LEDQ, DWQL, DUWQ, DWLQ, WDQE, 
DBUE, LGQU, WDUQ, DLQE, WBGQ, WLGQ, BDEU, LDUQ, BUDE, LQDU, 
DQWE, UWLQ, WGUQ, LDQU, LDQE, WLDQ, DLUQ, WBDQ, DWQU, WGBQ, 
UWQL, GWQU, DQLE, GWLQ, WQGL, GLUQ, DWBQ, GLQU, WLUQ. (A.17) 


Additionally, as shown in table 4, there are 176 newly forbidden 5-particle hierarchies, 
426 newly forbidden 6-particle hierarchies, 434 newly forbidden 7-particle hierarchies, and 
22 newly forbidden 8-particle hierarchies. In order to save space, we do not list these 
hierarchies explicitly here, however they can be easily obtained from the python code 
described in appendix B. 
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A.8 Model m = 8 


In this model, all 2-particle and 3-particle hierarchies are allowed. The forbidden 4-particle 
hierarchies are: 

BGEU, GBWQ, WUQL, GEBU, GWQL, WBGQ, WGBQ, WGLQ, 

BEGU, WGQL, GBEU, UWQL, BGUE, GWQB, WUGQ, GUWQ, 

UWLQ, WGQB, WGUQ, WLUQ, UQWL, WGQU, GWBQ, GWUQ, 

GBUE, GWQU, GWLQ, WULQ, GLWQ, WQGL, WLGQ, WQUL. (A.18) 

Additionally, as shown in table 4, there are 148 newly forbidden 5-particle hierarchies, 
809 newly forbidden 6-particle hierarchies, 398 newly forbidden 7-particle hierarchies, and 
54 newly forbidden 8-particle hierarchies. In order to save space, we do not list these 
hierarchies explicitly here, however they can be easily obtained from the python code 
described in appendix B. 

B Python package for studying hierarchies 

We have provided a Python [82] package in the supplemental material to the arXiv ver¬ 
sion of this paper. It consists of two required files, susyJiierarchyjnethods .py, and 
allowed_hierarchies .py and two optional files, f orbidden_hierarchies .py and 
newly_f orbiddenJiierarchies .py. These files should work with all versions of Python 
including and subsequent to Python 2.6. 

Information about which hierarchies are allowed or forbidden in various SUSY scenarios 
is contained in the files allowed_hierarchies.py and forbidden_hierarchies.py. Of 
course a hierarchy in a given SUSY scenario is either allowed or forbidden, so it is sufficient 
to have only, e.g., an enumeration of allowed hierarchies in various models. However, 
the generation of forbidden hierarchy information from allowed hierarchy information is 
relatively slow ~ 1 minute on a modern laptop, so we do include the, in principle redundant, 
files forbidden_hierarchies.py and newly.forbidden_hierarchies.py. 

B.l Data structures in Python 

To explain how hierarchy data is stored in the files allowedTiierarchies.py, 
f orbidden_hierarchies .py, and newly_forbidden_hierarchies .py, we must briefly re¬ 
view a few of the data structures present in Python. A list consists of elements, in a fixed 
order, between the square brackets ‘[’ and ‘[’. Elements are indexed by integers, e.g., the 
first element of niy_list is my_list [0]. A dictionary is like a list, but the elements are 
stored between curly brackets ‘{’ and and are indexed by arbitrary items called “keys”. 
An element of a dictionary is accessed by its key. For example if 

my_dict = {'example37, 32.3different'} 

then 

my_dict['example'] = 37 
my_dict['32.3'] = 'different' 
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Finally a tuple is like a list, except that tuples cannot be modified (and tuples are enclosed 
in parentheses). 

B.2 Hierarchy Data 

The file allowed_hierarchies .py contains the single command 

allowed = {2 : {('L', 'W'): [1, 3, 5, 2, 7, 4, 6, 8], ('Q','G'): \ 

[1, 3, 5, 2, 7, 4, 6, 8], .... 

The dictionary “allowed” has keys 2, 3, ..., 8; these keys represent the number of particles 
in a hierarchy. The elements corresponding to each of these keys are dictionaries. The keys 
for these dictionary are the hierarchies, which are tuples of capital letters representing the 
particles in the hierarchy. The entry corresponding to a hierarchy is a list (in no particular 
order) of the SUSY scenarios in which the hierarchy is allowed. The numbering of SUSY 
scenarios or “cases” is that which is used throughout this paper and defined in table 2. 
Note that only hierarchies which are allowed in some hierarchy have entries in “allowed”. 
Thus 

print(allowed[2][('U','Q')]) 
yields 

[3, 5, 7, 4, 6, 8] 

as the 2-particle hierarchy ('U' , 'Q') is allowed in models 3, 4, 5, 6, 7, and 8, while 

print(allowed[4] [('B', 'G', 'E', 'U')]) 

yields the error (in interactive session): 

Traceback (most recent call last): 

File "<stdin>", line 1, in <module> 

KeyError: ('B', 'G', 'E', 'U') 

as the 4-particle hierarchy ('B', 'G', 'E', 'U')is not allowed in any model. 

The file f orbidden_hierarchies .py uses the exact same format, but the dictionary is 
called “forbidden” and contains lists of the models in which hierarchies are forbidden, rather 
than allowed. Likewise the file newly_f orbidden_hierarchies .py uses a dictionary called 
' 'newly_forbidden' ', which tells us which hierarchies are newly forbidden (irreducible). 

B.3 Methods 

The Python expert can use the data in allowed Jiierarchies . py, forbidden Jiierarchies. py, 
and newly_f orbiddenJiierarchies .py for their own purposes. However, for those less fa¬ 
miliar with Python or for those who would benefit from an example of methods (functions) 
which use the hierarchy data, we have provided another module (python file) with meth¬ 
ods that process the data in allowed_hierarchies.py, forbidden_hierarchies.py, and 
newly_f orbidden Jiierarchies .py. These methods are 
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1. convert_to_standard_form(s) 

This method translates a string, s, to a tuple of capital letters representing the SUSY 
particles in a hierarchy. This is the format of the keys representing hierarchies in our 
data arrays. 

2. print.case(case) 

This method prints a string describing the SUSY scenario (“case”) labelled by the 
integer case, which must be between one and eight (inclusive). 

3. get_info_on_hierarchy(s) 

This method prints information listing and describing the SUSY scenarios (“cases”) 
in which the hierarchy s (s is a string, enclosed by single or double quotes) is allowed 
and in which it is forbidden. 

4. number_of_N_hierarchies_in_case(data_type,N,case) 

This method returns the number of N particle hierarchies which are either allowed or 
forbidden (depending on whether the value of datatype is “allowed” or “forbidden”) 
in SUSY scenario case. 

5. list_of_N_hierarchies_in_case(data.type,N, case) 

This method returns a list of N particle hierarchies which are either allowed or 
forbidden (depending on whether the value of datatype is “allowed” or “forbidden”) 
in SUSY scenario case. 

6. print_table_of.cases 0 

This method prints a table that describes the eight SUSY scenarios (“cases”) that 
we consider. 

7. print_table_ofjiumber_of.hierarchies (data.type) 

This method prints a table of all particle hierarchies which are either allowed or 
forbidden (depending on whether the value of datatype is “allowed” or “forbidden”) 
in all eight SUSY scenarios. 

B.4 Using the code 

The Python expert can use the methods and as part of their own code by importing 
whichever modules they need. However, the simplest way to use the methods provided is 
to use Python in interactive mode. To do this, one types 

python 

in a terminal. While it is discouraged in general (as it negates useful features of Python in 
regard to namespaces), the simplest next step is to type 

from susy_hierarchy_methods import * 

into the Python interpreter. Note that the files susyJiierarchyjnethods .py and 
allowed_hierarchies .py must be in the current working directory, or must be in one’s 
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PYTHONPATH. Doing this also imports “allowed” from allowedJiierarchies .py automati¬ 
cally. If forbiddenJiierarchies .py is in the working directory or the PYTHONPATH, then 
“forbidden” will be imported from this file. If not, “forbidden” will be generated automati¬ 
cally, though this can be somewhat time-consuming (~ 1 minute on a modern laptop). The 
situation regarding newly_forbidden_hierarchies .py is exactly analogous. One can then 
simply type the method one wishes to execute together with the appropriate arguments. 
Finally, we note that typing 

help(naine_of _niethod) 

will provide a brief description of the method in question. 
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